Fractional Fermions with Non-Abelian Statistics 
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We introduce a novel class of low-dimensional topological tight-binding models that allow for 
bound states that are fractionally charged fermions and exhibit non-Abelian braiding statistics. The 
proposed model consists of a double (single) ladder of spinless (spinful) fermions in the presence 
of magnetic fields. We study the system analytically in the continuum limit as well as numerically 
in the tight-binding representation. We find a topological phase transition with a topological gap 
that closes and reopens as a function of system parameters and chemical potential. The topological 
phase is of the type BDI and carries two degenerate mid-gap bound states that are localized at 
opposite ends of the ladders. We show numerically that these bound states are robust against a 
wide class of perturbations. 

PACS numbers: 71.10.Fd; 05.30.Pr; 71.10.Pm 



Introduction. Topological properties of condensed 
matter systems have attracted considerable attention in 
recent years. In particular, Majorana fermions [l|, be- 
ing their own antiparticles, are expected to occur in a 
number of systems, e.g. fractional quantum Hall systems 
EEC! , topological insulators optical lattices pill], p- 
wave superconductors J9| , nanowires with strong Rashba 
spin orbit interaction [lCHla ] , and carbon-based systems 



16l-ll8l|. Another class of topological systems is given by 



bound states of Jackiw-Rebbi type [l9[ containing frac- 
tional charge e/2 [22 - 27 1 . Such exotic quantum states are 
interesting in their own right, and due to their special ro- 
bustness against many forms of perturbations they offer 
the possibility for applications in quantum computations, 
especially when they exhibit non-Abelian statistics such 
as Majorana fermions 2fj, 21 1. 




In this Letter we introduce a surprisingly simple class 
of models supporting a topological phase with bound 
states that possess not only fractional charge but also ex- 
hibit non-Abelian statistics under braiding. These bound 
states behave in many ways similar to the well-studied 
Majorana fermions in superconducting-semiconducting 
nanowires [20[ , but in contrast to them they are complex 
fermions, and quite surprisingly, emerge in the absence 
of superconductivity and without BCS-likc pairing. 

The two non-interacting tight-binding models we pro- 
pose consist of a double ladder containing spinless parti- 
cles in a uniform magnetic field and a single ladder con- 
taining spinful particles in the presence of both uniform 
and spatially periodic magnetic fields. We find a topo- 
logical phase transition in these systems when varying 
system parameters or the chemical potential, with a char- 
acteristic closing and re-opening of a topological gap. In- 
side the topological phase we find two degenerate bound 
states, one localized at the right and one at the left end of 
the system. These bound states are fractionally charged 
fermions and are shown to exhibit non-Abelian braiding 
statistics of the Ising type. We study the systems analyt- 
ically in a continuum approach, finding explicit solutions 



FIG. 1. Double-ladder tight-binding model, consisting of a 
lower (a — —1) and an upper (a = 1) ladder, each lying 
in the xy-plane and held at chemical potentials /i CT . Here, 
t x (red links) is the intra-chain, t y (green links) the intra- 
ladder, and t z (blue links) the inter-ladder hopping amplitude, 
and the corresponding lattice constants. A uniform 

magnetic field B = Bi + B2 is applied in the j/z-plane. The 
associated magnetic flux results in phases <f>i and 4>2 in the 
hopping amplitudes between different chains, see Eqs. (|2I3[) . 



for the bound states, and confirm these findings by inde- 
pendent numerics of the underlying tight-binding model. 
We further test the stability of these states numerically 
against a wide class of perturbations and show that the 
bound states are robust against most of them, except of 
local charge fluctuations, against which they are partly 
protected by charge neutrality. 

Tight-binding model. We consider a double-ladder sys- 
tem consisting of four coupled chains aligned along x- 
direction, see Fig. [TJ Two upper (lower) chains form 
the upper (lower) ladder. Each chain is labeled by two 
indices r and a, where r = ±1 refers to the left /right 
chains, and a = ±1 refers to the upper /lower ladders. 
The tight-binding Hamiltonian of a (r, er)-chain reads 



cr,n+l i",(T,rt 



h.c. 



(M L T,IJ,II, 



(1) 
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where c^o-.n (c\. an ) is the annihilation (creation) opera- 
tor on site n of the (t, er)-chain. The sum runs over N 
sites composing the chain. Here, t x is the hopping matrix 
element in x-direction, and \i a is the chemical potential 
of the (r, <r)-chain. 

The intra-ladder coupling is given by 

HI = t v Yy n ^ c \^x,«,n + h.c), (2) 

n 

where the phase </>i accompanies the hopping matrix el- 
ement t y in ^-direction. This phase arises from the mag- 
netic flux through the unit cell produced by a magnetic 
field Bi in z-direction. The inter-ladder coupling be- 
tween two left or two right chains is given by 

H* = t t ^(e in ^4 il)n c TiIin + h.c). (3) 

n 

Similarly, the phase 4>2 arises from the flux through 
the unit cell produced by a magnetic field B 2 in y- 
direction. We note that in practice only one total field, 
B = Bi + B2, needs to be applied in the yz-planc, 
sec Fig. [T] The total tight-binding Hamiltonian for the 
double-ladder model is given by 

H = Y}H x r , a + m + H*)l% (4) 

Now we focus on a particular case of the above model. 
First, we fix the chemical potentials on the upper and 
lower ladders to be of opposite signs, i±\ = —\x\. Second, 
we assume the upper ladder to be at quarter-filling, i.e. 
Ill = s/2t x . The magnetic fields are chosen such that 
<±>i = 7r/2 and (f>2 — 7r, or in terms of field strengths, 
B% = §o/Aa x a y and B2 = $o/2a x a z , where $0 is the 
flux quantum, and a XtVlZ are the corresponding lattice 
constants. Assuming that t x 3> t y ,t z , we treat from now 
on inter-chain hoppings as small perturbations. 

Continuum model. The most convenient way to an- 
alyze the tight-binding Hamiltonian in Eq. ((4]) is to 
go to the continuum limit. For this we first derive 
the spectrum via Fourier transformation along the x- 
axis for each (t, cr)-chain, c TjCTi „ = -±= J2k elkna!C c r,<7,k, 
where c r (T jt is the annihilation operator of the electron 
with momentum fc. The Hamiltonian a has the well- 
known spectrum, €k,a,n = AV + 2t x cos(ka x ). At quarter- 
filling = ±V2t x ) the Fermi momenta are given by 
kp.i = ir/4a x and k F \ = 3n/4a x , see Fig. [2] We em- 
phasize here that the system is charge neutral, and the 
shifted chemical potentials just redistribute electrons be- 
tween chains. 

Next, we linearize the spectrum around the Fermi mo- 
menta by expressing the annihilation operators ty(x) that 
act on the states close to the Fermi level in terms of slowly 
varying right (R T ^ a ) and left {L T , a ) moving fields as 

*(a0 = [e lkF -° x RrAx) + e- tk ^ x L T ^(x)] . (5) 

T,<7 
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FIG. 2. The spectrum of the upper (orange) and lower (red) 
chains in the first Brillouin zone. The chemical potentials 
are chosen such that the system is at the charge-neutrality 
point. The filled (empty) states are indicated by dark-colored 
(light-colored) lines. The two Fermi wave vectors are given by 
kF.i = ir/4:a x and k F j = 3n/iax- The intra-ladder hopping 
t y (green dashed line) and the inter-ladder hopping t z (blue 
dotted lines) lead to opening of gaps at the Fermi level (e = 0) . 

The kinetic part of the Hamiltonian corresponding to 
J2 T a a [see Eq. (JTJ] is rewritten as 

H x = -ihv F J2 I da (Rl a d x R T , a - l..,J),l.:. : ) , (6) 

where we dropped the fast oscillatory terms, and where 
vp = \/2t x a x /h is the Fermi velocity. The intra-chain 
couplings, given by H* and H y [see Eqs. and (J3J], 
lead to mixing between i? Ti(T and £ Ti<T belonging to dif- 
ferent chains. The intra-ladder hopping yields 

H y = J dxt y (R\ tl Li tl + R\ tl L lt i + h.c)) , (7) 

while the inter-ladder hopping yields 

H z = J dxt z (R\ tl L^ + R} u L u 

■ />';. ,/.,.: I /',../•:. 1 • /'•<••)• (8) 

Next, we introduce a new basis <f> = 

to rewrite 

the total Hamiltonian H = H x + H y + H z as 
H = J dx (fy (xyhL<f){x) in terms of the Hamiltonian 
density W, 

H = hv F krj 3 + t z ai-qi + y (n»7i - ^0^2) + £m°3, (9) 

where the Pauli matrix rji acts on the right- and left- 
mover subspace, and the Pauli matrices o~i and act 
on the chain subspaces. The momentum operator is de- 
fined as k = —ihd x , with eigenvalue k counted henceforth 
from the corresponding Fermi points. Here, we assume 
again that the chemical potentials of the upper and lower 
chains are opposite in sign, however, small deviations 
from quarter-filling, 5fi = 5fXi = —Sm, are taken into 



3 



account. In addition, wc neglect any constant shifts of 
the spectrum. 

The Hamiltonian density % allows us to determine the 
topological class of the system [28j |. The system is in- 
variant under the time-reversal operation T, defined by 
UlU*(-k)U T = H(k). Indeed, U T = mn satisfies this 
relation. Similarly, the charge-conjugation symmetry op- 
eration C, defined by U^H*(-k)U c = -U{k), can be sat- 
isfied by Uc = T\0\f}%- Thus, the system belongs to the 
topological class BDI J^. The one-dimensional systems 
of this class are allowed to have an arbitrary number of 
bound states inside the energy bulk gap [28[. To deter- 
mine if bound states are present in the system for some 
given set of parameters, we follow the method developed 



in Rcfs. 27 andk2 



The eight spectrum branches of % are given by 

(ef 2 ) 2 = (hv F k±Sfi) 2 +t 2 , 
(e± 4 ) 2 = (hv F k) 2 +i 2 z +t 2 + d> 2 



(10) 



8fi 2 [(hv F k) 2 + t 2 }. (11) 



The system is gapless only for one particular set of pa- 



rameters, ti 



, 8fi . Otherwise, the spectrum is 
gapped, and there is a possibility for the existence of 
bound states inside this gap. Further, we are interested 
in states exactly in the middle of the gap, i.e. at zero en- 
ergy. In addition, wc focus on semi-infinite chains, so the 
boundary conditions are imposed only on the left (right) 
end. This implies that the chain length L is much larger 
than the localization length £ of the bound states we find. 

In order to address the existence of bound states, we 
first find four fundamental decaying solutions of the sys- 
tem of linear differential equations, following from the 
Schrodingcr equation associated with H [see Eq. ©J. 
Second, the dimension of the null space of the corre- 
sponding Wronskian leads us to a topological criterion 
[27L [29l | that separates a topological phase (with bound 
states) from a trivial phase (without bound states). We 
find that bound states exist provided the following topo- 
logical criterion is satisfied, 



t\ > t 2 z +5fi 2 . 



(12) 



Working in the operator basis ^(x) = 
(*i,1i^Hi^IIi*i i)i where ^ T)(T (x) is the anni- 
hilation operator on the (r, er)-chain, we find the 
wavefunction of the state localized at the left end of the 
double-ladder system explicitly, 



ib L {x) = 



/*(*) 
V -(-l) n /(x) / 



(13) 



and the wavefunction of the state localized at the right 
end, ipn( x ) = ^^{F — x), where n labels the site, x = na x 



N. 
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FIG. 3. (a) The relevant part of the spectrum for a double- 
ladder system of length L/a x — 401 found by numerical diag- 
onalization of the tight-binding Hamiltonian H [see Eq. Q]. 
The energy eat corresponds to the Nth. energy level. The pa- 
rameters t~/t x = 0.05, ty/tx = 0.1, and S/j, — are chosen to 
satisfy the topological criterion Eq. (|12|) . In the middle of the 
gap, e = 0, there are two degenerate bound states (red dots), 
(b) The probability density I^t.^I 2 of the left state in the 
(r, <r)-chain. The density patterns are the same for all chains, 
in full agreement with the continuum solution Eq. (|13[) . 



[3fJ| . Here, we have introduced the notations, 

(14) 



The localization length of the bound states is given 
by £ = max{£i, £2}, where hv F ^i = t z , and hv F £,2 = 
yjt 2 — (5/z 2 — t z . Close to the phase transition point the 
localization length is determined by £ 2 , whereas deep in- 
side the topological phase by £1. If £ is comparable with 
L, the two bound states localized at opposite ends over- 
lap. As a result, the energy levels are split from zero 
energy, and the corresponding wavefunctions are given 
by the symmetric and antisymmetric combinations, 

^s/a(x) = 1pL(XL) ± ilpR(xR). (15) 

We note that these right or left localized bound states 
are fractionally charged fermions of charge e/2 22h24|. 

The results obtained above in the continuum model are 
in good agreement with the numerical results obtained by 
direct diagonalization of the tight-binding Hamiltonian 
H given in Eq. (gj, see Fig. [3l 

Non-Abelian statistics. In the absence of overlap be- 
tween the bound states ?/>l(x) and i/jr(x), the zero-energy 
level is two-fold degenerate, so it can potentially be used 
for braiding and ultimately for topological quantum com- 
putation. Taking into account that the system is charge- 
neutral, we focus on the states that have equal density 
distributions on both ends, 



i>±{x)=^ L {x)±e ia i> R {x), 



(16) 



where a is an arbitrary phase (see also below). The cre- 
ation operator F± corresponding to the state tp±(x) is 
given by F± — J dx ip±(x) ■ * = -^{Fl ± e %a F R ), where 
Fl (Fr) is the left (right) endstate annihilation operator. 



During the braiding process, which corresponds to ex- 
changing the right and left end states, the parity, defined 
by P± = 1 — F±F±, should be conserved [21(. As a result, 
the old states transform into new ones as F± — > ±e~ la F± 
under braiding, or in terms of the left/right states, 



F L -> Fr, F r 



e~ 2la F L . 



(17) 



The corresponding unitary operator Ulr that imple- 
ments this braiding rule is found to be Ulr = 
e -ia(n + +n_)+i7rn_ ) whcre n± = fIf±. Indeed, it is 

easy to show that U\ jR F±Ulr = ±e~ ta F±, and thus 
U f LR F L U LR = F R , and Ul R F R U LR = e~ 2la F L . In terms 
of the Fr/l operators we have 

U LR = e -i(a—ir/2)(n R +n L )-iir(e lc "FlF R +h.c.)/2 r^g\ 



where 



n L/R - F L /R F L/R- 
work of such double ladders similar to the one proposed 
for Majorana fermions [3l| . The braiding operations can 
be performed by exchanging bound states localized at 
different double ladders. Using Ulr given explicitly by 
Eq. (fT8|) , one can show that two braiding operations do 
not commute, [Uij, Ujk] ^ 0, where the labels k de- 
note three states. Thus, we see that our bound states 
obey non-Abelian braiding statistics. In particular, it 
is interesting to consider the case of a = ±7r/2, since 
these states, ^ s / a {x) [see Eq. (|T5|) ]. can be easily pre- 
pared by lifting the degeneracy temporarily (via tuning 
the wave function overlap), resulting in filling either the 
symmetric or the antisymmetric energy levels after some 
dephasing time. Moreover, we can determine the phase a 
[see Eq. ([To]) ] by projective measurements on the states 
V's/o- F° r example, the probability to measure the sym- 
metric state is given by | (ip s \?p + ) | 2 = sin 2 (a/2 + 7r/4). 
All this taken together opens up the possibility to use 
these bound states for topological quantum computation 



Next, let us assume a net- 



along the lines proposed for Majorana fermions [31 [ . 

Single-ladder model. An alternative representation 
of the double-ladder model with spinless particles is a 
single-ladder model with spinful particles, see Fig. [4] In 
this model the a = — 1 (lower) and the a = 1 (upper) lad- 
ders are identified with the spin-down and spin-up com- 
ponents, resp. A uniform magnetic field Bi is applied 
perpendicular to the ladder, leading to both orbital and 
spin effects. The Zeeman energy acts as a chemical po- 
tential with opposite signs for opposite spin directions, 
fi a = agfiRBi, where g is the ^-factor, and fiR is the 
Bohr magneton. A spatially periodic field §2 of period 
2a x couples opposite spins, so the effective hopping is 
given by t z = g^B^2- Such a periodic field can be pro- 
duced by nanomagnets [32| or, equivalently, by Rashba 
spin orbit interaction and a uniform magnetic field (33| . 

Stability against perturbations. We next address the 
stability of the topological phase and the bound states 
against local perturbations. In general, we find numer- 
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FIG. 4. Single-ladder model. A uniform magnetic field Bi 
is applied perpendicular to the ladder. A spatially periodic 
magnetic field B2 with period 2a x is applied in the xy-pl&ne. 



ically that local fluctuations that preserve the symme- 
try between the upper and lower ladders are not harm- 
ful. This includes correlated fluctuations of chemical po- 
tentials (A/zi = — A/ij), magnetic fluxes, and hopping 
matrix elements. If the symmetry is not preserved, the 
perturbations act like a level-detuning, and the bound 
states separate in energy independent of wave function 
overlap but proportional to their occupation probabil- 
ity at the site of fluctuation. We emphasize that the 
single-ladder model is protected against such symmetry- 
breaking terms, except against local chemical potential 
fluctuations that break charge neutrality, i.e. A/i! = 
Afij. However, we note that chains without charge im- 
purities are stable against such fluctuations, as local dif- 
ferences in /i would lead to charge redistribution, restor- 
ing a uniform chemical potential in the chain. Another 
problem can arise from flux fluctuations. They can de- 
crease the Fourier components at kF,a of the backscat- 
terig terms and thereby reduce the gaps. As a result, the 
system can move out of the topological phase. However, 
these flux fluctuations become irrelevant deep inside the 
topological phase. 

Conclusions. We have uncovered model systems of 
striking simplicity that allow for a topological phase with 
degenerate bound states that are fractionally charged 
and obey non-Abelian braiding statistics. We have 
shown that these exotic quantum states are rather robust 
against a large class of perturbations. Quite remark- 
ably, our models demonstrate that non-Abelian states 
can exist in single-particle systems, without any correla- 
tions and in the complete absence of superconductivity or 
BCS-likc pairing. This should open the path for novel im- 
plementations of topological matter in realistic systems. 
One promising candidate system that suggests itself for 
implementations of such tight-binding ladders are opti- 
cal lattices [Hj], because they allow for a high degree of 
control and, in particular, possess the charge stability of 
the type invoked here. 

This work is supported by the Swiss NSF, NCCR 
Nanoscience, and NCCR QSIT. 
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